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Catalysis of dynamical symmetry breaking by a constant magnetic field in (3 + 1) dimensions is 
considered. We use the three flavour Nambu - Jona-Lasinio type model with 't Hooft and eight-quark 
interaction terms. It is shown that the multi-quark interactions introduce new additional features 
to this phenomenon: (a) the local minimum of the effective potential catalyzed by the constant 
magnetic field is smoothed out with increasing strength of the field at the characteristic scale 
H ~ 10 19 G, (b) the multi-quark forces generate independently another local minimum associated 
with a larger dynamical fermion mass. This state may exist even for multi-quark interactions with a 
subcritical set of couplings, and is globally stable with respect to a further increase of the magnetic 
field. 



PACS numbers: 11.30.Rd, 11.30.Qc 

It has been shown in a series of papers [l|, 0, 9 that in 
(2 + 1) and (3+1) dimensions a constant magnetic field 
H 7^ catalyzes the dynamical symmetry breaking lead- 
ing to a fermion mass even at the weakest attractive four- 
fermion interaction between particles, and the symmetry 
is not restored at any arbitrarily large H . Soon thereafter 
it became also clear 0, [H, H| that the zero-energy surface 
of the lowest Landau level (LLL) plays a crucial role in 
the dynamics of such fermion pairing. It has been found 
that the dynamics of the fermion pairing in the homo- 
geneous magnetic field is essentially (1 + l)-dimensional, 
and a deep analogy of this phenomenon with the dynam- 
ics of electron pairing in BCS Q has been stressed. The 
generated fermion mass, Md yn , turned out to be much 
smaller than the Landau gap ~ y|eiT|. 

The existence of a zero-energy surface in the spectrum 
of a Dirac particle is ensured for any homogeneous mag- 
netic field with a fixed direction by a quantum mechan- 
ical supersymmetry of the corresponding second-order 
Dirac Hamiltonian Q. This aspect of the phenomenon 
appears to be a quite exceptional situation and indicates 
that the dynamical generation of mass is not so univer- 
sal as one would expect by extrapolating the results ob- 
tained for homogeneous or unidirectional Q magnetic 
field profiles. For instance, it has been demonstrated 
by Ragazzon that the Nambu - Jona-Lasinio (NJL) 
model [lll | minimally coupled to a background magnetic 
field with variable direction does not possess a massive 
phase until the coupling constant exceeds some critical 
value. Obviously, in this case one faces the conventional 
scenario of dynamical chiral symmetry breaking, where 
the magnetic field does not play an essential role. 

Conversely, having in mind that homogeneous mag- 
netic fields can act as strong catalysts of chiral symme- 
try breaking, one might ask what is the effect caused by 
the strong interaction, when higher order multi-fermion 
interactions are present. These extensions of the NJL 



model are well-known [T^, [H, [H| j for instance, the four- 
quark U(3)l x U(3)r chiral symmetric Lagrangian to- 
gether with the U(1)a breaking 't Hooft six-quark in- 
teractions has been extensively studied at the mean-field 
leve l fl5l [H, [I?], 13- Recently it has been also shown 
[H^ , |20|] that the eight-quark interactions are of vital im- 
portance to stabilize the multi-quark vacuum. 

The additional multi-quark forces can affect the result 
which is obtained when only four-fermion interactions are 
considered. We argue, in particular, that the 't Hooft 
and eight-quark interactions can modify the theory in 
such a way that the local minimum, catalyzed by the 
constant magnetic field, is smoothed out by increasing 
the strength of the field. This is an alternative regime to 
the known one in which the strong magnetic field cannot 
wash out the condensate from the LLL. For the first sce- 
nario to become possible it is sufficient that the couplings 
of multi-quark interactions are chosen such that the sys- 
tem displays more than one solution of the gap equation 
at H = 0. Actually, the above condition is not a require- 
ment. Even if the gap equation has only one nontrivial 
solution at small H , an increase in the magnetic field can 
induce the formation of a second minimum. 

The multi-quark dynamics of the extended NJL model 
is described by the Lagrangian density 

£ G ff = qii^D^ - m)q + C iq + C 6q + C Sq , (1) 

where the gauge covariant derivative is equal to 
Dfi = (9 M + iQA^ with the external electromagnetic field 
and quark charges Q = e ■ diag(2/3, —1/3, —1/3). It 
is assumed that quark fields have colour (N c = 3) and 
flavour (Nj = 3) indices. The current quark mass, m, is 
a diagonal matrix with elements diag(m u , rhd, rh s ), which 
explicitly breaks the global chiral SUl(3) x SUr(3) sym- 
metry of the Lagrangian. We shall neglect this effect in 
the following assuming that rh = 0. 

The multi-quark interactions (in the scalar and pseu- 
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field H. As has been calculated by Schwinger a long time 
ago [13 



where 



V s = Vs(mi,\QiH\), 



(10) 



The U(3) flavour matrices A a , a = 0, 1, ... ,8, are nor- 
malized such that tr(A a Ab) = 25 a b- The matrices Pl,r — 
(1 =p 75)/2 are chiral projectors and the determinant is 
over flavour indices, which are suppressed here. The de- 
terminantal interaction breaks explicitly the axial U(1)a 
symmetry [2l| and Zweig's rule. The eight-quark spin 
zero interactions are given by 

£& g = 8 5i [{q~iPRqm)(q m PLqi)] , (5) 
A ? = ^92(qiPRq m ){q m PLqj){qjPRqk){qkPLqi)- (6) 

G, k, gi, g 2 are dimensionful coupling constants: [G] = 
M~ 2 , [k] = M~ 5 , [gi] = [g 2 ] = M~ 8 in units h=c=l. 

We proceed by calculating the effective potential of 
the theory, V(m tt) m<j, m s ), in a constant magnetic field: 
A x = —Hy, A y = A z = (Landau gauge). The argu- 
ments, rrii, are simply real parameters; they are not to be 
identified with the masses of any presumed one-particle 
states. Instead, we shall use the capital letter Mj for the 
point where V takes its local minimum, which specifies 
the masses of constituent quark fields. 

The potential is built of the following two terms 



V(m u , m d , m s ) = V st + V s . 



(7) 



The first contribution results from the many-fermion ver- 
tices of Lagrangian £ e ff > after reducing them to a bilin- 
ear form with help of bosonic auxiliary fields, and subse- 
quent integration over these fields, using the stationary 
phase approximation (SPA) method. The specific details 
of these calculations and the result are given in a recent 
work [lj|. We obtain 

V st = i (±Gh 2 + nh u h d h s + ^ (h 2 ) 2 + 3g 2 h^j , 

(8) 

where h 2 = V, , iti Jr. and /i 4 = E i=u , d ,A-- The 
functions hi depend on the coupling constants G, k, g±, g 2 
and on the independent variables A, = m$ — rhi. To 
find this dependence one should solve the system of cubic 
equations 



Ghi 



h Ir + y ^ = 0. (9) 



In some parameter range the system has only one set of 
real solutions, and this guarantees the vacuum state of 
the theory to be stable [19| ■ 

The second term on the r.h.s. of eq. ([7]) derives from 
the integration over the quark bilinears in the functional 
integral of the theory in presence of a constant magnetic 
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' p{s,A 2 )\QH\coth(s\QH\) + const. 



Here the cutoff A has been introduced by subtracting 
off suitable counterterms to regularize the integral at the 



lower limit, i.e., p(s,A ) = 1 — (1 



sA 2 



For the 



fcrmion tadpole this works as the four-momentum co- 
variant cutoff in the euclidean space: p 2 +p 2 < A 2 . The 
unessential constant is chosen such that Vs(0, \QH\) = 0. 
As a result we obtain 
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The quantity C(y , x) denotes the generalized Riemann 
zeta function [231 ] . 

We shall now illustrate the procedure which will be 
employed in the following, by considering first the simple 
SU(3) flavour limit for the situation in which m = and 
n = gi = g 2 = 0. For the purpose of illustration, we 
ignore the charge difference of u and d, s quarks in the 
remaining. The averaged common charge \Q\ — |4e/9| 
will be used. In this case one obtains the potential 
Vim) = N f (m 2 /4:G + V s (m, \QH\)). One sees that the 
gap equation, dy(m)/dm = 0, has always a trivial solu- 
tion, tyl — 0. The nontrivial solution is contained in the 
equation 

27t2 tl 2 a ™x , M 2 + m 2 \ \QH\ 

2 In \-'-'- | - / , (13) 
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where ip{ x ) = dlnT(a;)/da; is the Euler dilogarithmic 
function. This equation has a solution at all G > 0, 
if H ^ 0. Fig. [T] illustrates this important result for 
GA 2 = 3. One sees that in absence of the magnetic 
field the system is in the subcritical regime of dynami- 
cal symmetry breaking. The introduction of a constant 
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asymptotic behaviour is 



FIG. 1: The l.h.s. (straight short-dashed line) and the r.h.s. 
of eq. (|13[) as functions of m/A for four different values of 
the magnetic field strength H: full curves (top to bottom) 
correspond to |Q-ff|A~ 2 = 0.5; 0.3; 0.1, and the dashed curve 
to H — 0. Box insert: close-up of region around origin with 
solid lines for IQJ^IA -2 = 0.2; 0.15; 0.1 (top to bottom). 



field, however small it might be, changes radically the 
dynamical symmetry breaking pattern, due to the singu- 
lar behaviour of the r.h.s. of cq. (fT3")) close to the origin: 
the right and left hand sides will always intersect and 
the value of m where this happens is a minimum of the 
effective potential. 

Let's look at these different regimes more closely. As 
H — > 0, we recover the well-known NJL model gap equa- 
tion 



1 = 



GN C 
2vr 2 



Jo(m 2 ) 



\QH\ 2 A 2 



3m 2 (A 2 + m 2 ) 



(14) 



where Jo(m 2 ) = A 2 - m 2 In (l + A 2 /m 2 ) . Eq. (fT5J) at 
H = admits a nontrivial solution only if r > 1, where 
t = GA 2 N c /2ir 2 . This determines the critical value 
G cr = 2^ 2 /A 2 7V c . 

At m 2 /A 2 < 1 the r.h.s. of eq. (US]) is 
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Here the function v(£) of the argument £ = A 2 /2\QH\, 
is given by 

w(0 = ^[l-21nr(e + l)]+V(0- (16) 

This is a monotonically increasing function on the inter- 
val < £ < oo; v(£) = at the point £ ~ 1.12; the 
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L--ln(27rO, (£ - oo), (17) 

«(0~-7-^. (£-0), (18) 

where 7 ~ 0.577 is the Euler's constant. 

In the approximation considered one finds the solution 
of eq. (fH 



Mdyn = Aexp 



(19) 



To discuss the physical content of this result, we re- 
call that the energy spectrum of relativistic fcrmions in 
a constant magnetic field H contains Landau levels 



E n { Pz ) = ±v/m 2 + 2\QH\i 



Pz 



n = 0,1,2, 



(20) 



with p z denoting the projection of the 3-momentum on 
the z-axis, i.e., along the magnetic field. If the fermion 
mass 771 goes to zero, as in the present case, there is no 
energy gap between the vacuum and the LLL. Thus the 
integer part of £ + 1 gives approximately the number of 
Landau levels taken into account. 

The first term in eq. (|15p has a clearly defined two- 
dimensional origin, given the logarithmic dependence on 
the cutoff in the corresponding gap equation 



^T |0H|ln A 2 " 



(21) 



and, therefore, in the condensate (compare with eq. 
(|14[) ). Such behaviour is associated with the (1 + 1)- 
dimensional dynamics of the fermion pairing on the en- 
ergy surface Eq = of the LLL Q. As long as this 
term dominates over the second term, v(£) in (fTS)) . one 
concludes that the condensate is mainly located on the 
LLL. Actually this condition is fulfilled nearly every- 
where at t < 1. Indeed, this is obvious for £ = 1, since 
v(l) = 1/2 + if>(l) = 1/2 - 7 ~ -0.08 is small compared 
with 1/t. For £ < 1 we come to the same conclusion af- 
ter considering the asymptotics of the second term (|18[) . 
The other formula, (|17[) . can be used to show that the 
above statement is also true for £ > 1, except near the 
critical region t — * 1 — 0, where v(£) dominates; then the 
condensate spreads over many Landau levels. 

In this special case it is possible to find an analytical 
solution. Indeed, using (fl~7|) in cq. (fT5|) we obtain 
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To progress further we suppose that the two following 
small variables are of the same order 



m 2 (\QH\y 

A 2 ~ I A 2 J ~ £ ' 



(23) 
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Then it follows immediately that the term with the loga- 
rithm is of order y/eln y/e and goes to zero, when e — > 0. 
Thus, the gap equation 



1 \QH\ fnm 2 \ 
1 -r = ^{\QH\) 



0(e) 



(24) 



is valid only in the region near the critical value r — ► 1—0. 
The closer r to 1, the smaller is e; Landau levels ap- 
proach a continuun distribution, and a condensate occu- 
pies many levels. The physical reason for the changes 
found in the behaviour of the condensate is the strength 
of the four-fermion interaction which becomes essentially 
important here. The corresponding solution is 



\QH\ 



exp 



A 2 fl _ 

]QH\ \t 



(25) 



Note that the near-critical regime found here, differs 
from the result of Ref. [5j, being driven by a quadratic 
dependence on the cutoff, eq. ([24]) ; this is tantamount of 
having a (3+l)-dimensional dynamics of fermion pairing. 

Let us return now again to the three flavour case with 
k, gi, <?2 7^ 0. In the simplest case with the octet flavour 
symmetry, where current quarks have equal masses m u = 
rhd = m s , which we set again zero, the system Q reduces 
to a cubic equation with respect to h = h u = hd = h s 



K ,2 4G 



4to 
3A 







(26) 



with A = g\ + (2/3)52- This cubic equation has one real 
root, if (see [20] for more details) 



G f k \ 2 
A > V24A/ 



(27) 



Assuming that the couplings fulfill condition (|27|) , we find 
a single valued function h(m) from eq. (|26[) . 

Considering that most of the investigations have been 
using multi-quark Lagrangians without the stabilizing 
eight-quark interactions we make a short digression to 
discuss the case with A = (see details in e.g. (l9l. |24|). 
In this case eq. ([26)) is quadratic with a regular and a 
singular solution as k — ► 0, i.e., 



= - 



8G 



urn 

Hp I 



(28) 



In SPA this leads to an unstable effective potential. In 
the commonly used mean field approximation, which 
discards the singular solution, the effective potential is 
metastable and the region 4G 2 < Km leads to complex 
values for h. This translates to a restriction for the ad- 
missible values that the l.h.s. of eq. (|29|) can assume, 
shown as dash-dotted line in fig. [2J discussed below. Fur- 
thermore we find that qualitatively the symmetry break- 
ing pattern is the same as in presence of the eight-quark 
interactions, but the occurrence of two minima requires 
higher values of \n\, as compared to the case with A 7^ 
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FIG. 2: The l.h.s. and the r.h.s. of eq. l[29)l as functions of 
m/A with /(m 2 ; A, \QH\) shown for |Q#|A- 2 = 0; 0.1; 0.5. 
The l.h.s. with A / (short-dashed line) is shifted with re- 
spect to the ordinate axis due to the 't Hooft interaction, the 
range of values for the A = curve (dash-dotted) is limited 
(see discussion in text) . The value of GA 2 = 3 is the same as 
in fig. [T] Here the six- and eight-quark couplings are taken to 
be kA 5 = -10 3 and AA 8 = 3.67 • 10 3 (or A = 0), respectively. 



(~ 2.2k for the parameter set of fig. [2]). From now on 
we shall consider only the case which fulfills the stability 
requirement, provided by eq. ([2T| . 

The nontrivial solutions for the mass of the fermion 
field in a constant magnetic field are determined by the 
equation 



2ir 2 h(m) 
A 2 N c m 



f(m 2 ;A,\QH\). 



(29) 



Comparing this result with eq. (fl"3)) . one sees that only 
the l.h.s. is changed. The six- and eight-quark interac- 
tions have modified it in such a way that now we get 
a function h(m)/m instead of the former constant term 
involving only the coupling of four-quark interactions, 
— 1/G. The l.h.s. of eq. ([2"9"| . abbreviated by u in fig. 
2, has now a bell-shaped form (short-dashed line), to be 
compared with the horizontal line of fig. 1. Note, that 
h(m)/m = —1/G + G(m), i.e., the bell-shaped curve 
crosses the ordinate axis at the same point as the former 
straight line (for the same value of GA 2 ). The r.h.s. of 
the eq. ([2"9"| is again represented by the long-dashed curve 
(H = case) and by the full lines for the finite H cases: 
\QH\A~ 2 = 0.1, 0.5. As mentioned before, these are not 
altered by the couplings G, k, A. The intersection points 
of the l.h.s. with the r.h.s. curves yield the nontrivial 
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solutions of the gap equation: one sees that either one or 
three solutions can be found for m > 0. If eq. (f2"9")l has 
no solutions at H = we say that the set of couplings 
G, k, A are subcritical. It is said to be overcritical in the 
opposite case. Note that the overcritical set may contain 
G < G cr . 

The trivial solution, m — 0, corresponds to the point 
where the potential, V(m), reaches its local maximum. 
Indeed, the second derivative 

hm^.lim^ln H^.-oo (30) 
nwo dm 2 m^a 2tt 2 2\QH\ v ' 

is negative here. This is the general mathematical rea- 
son for the phenomenon known as magnetic catalysis of 
dynamical flavour symmetry breaking. The logarithmi- 
cally divergent negative result ensures that this phase 
transition always takes place, if if ^ 0. This does not 
depend on the details related with the multi-quark dy- 
namics, i.e., the result is true even for free fermions in a 
constant magnetic field. 

What is really sensitive to the multi-quark dynamics 
is the local minima structure of the theory. Let us re- 
call that in the theory with just four-fcrmion interactions 
the effective potential has only one minimum at m > 0, 
and this property does not depend on the strength of 
the field if. We have demonstrated this in fig. [TJ In 
the theory with four-, six-, and eight-quark interactions 
one can find either one or two local minima at m > 0. 
The result depends on the strength of the magnetic field 
if, and couplings G, k, A. We illustrate these two cases 
in fig. [5J Namely, the upper full curve (r.h.s. of eq. 
(|29|) for |Qff |A~ 2 = 0.5) has only one intersection point 
with the bell-shaped curve u (l.h.s. of eq. (f2"9")l for 
GA 2 = 3, kA 5 = -1000, AA 8 = 3670). This point cor- 
responds to a single vacuum state of the theory. The 
other full curve (r.h.s. of eq. ([29]) for |C;if|A~ 2 = 0.1) 
has three intersections with the same curve u. These in- 
tersections, successively, correspond to a local minimum, 
a local maximum and a further local minimum of the 
potential. 

It is interesting to note that the first minimum cat- 
alyzed by a constant magnetic field (that is, a slowly 
varying field) is then smoothed out with increasing if. 
It ceases to exist at some critical value of |Qif |A -2 , from 
which on only the large M^ yn solution survives. This is 
shown in fig. [3J for the parameter set of fig. O This 
process is accompanied by a sharp increase in depth of 
the effective potential at the second minimum, especially 
if we had at the beginning the opposite ordering, i.e., 
V(Mi) < V{M2j- The reasons for such a synchronized 
behaviour are the following two facts. The first one is 
eq. (|30l) . which teaches us that the only way to wash 
out the first minimum is by lowering the barrier between 
this state and the second minimum. The other fact is 
the observation that the second minimum is unremov- 
able, because the asymptotic behaviour of the functions 
in eq. (|2T)|) is such that the l.h.s. dominates over the 
r.h.s. at large m/A. This can also be understood from 




0.05 0.1 0.15 0.2 



\QH\/A 2 

FIG. 3: The dimensionless dynamical mass Mdyn/A as a func- 
tion of the dimensionless magnetic field |Qif|/A 2 . The full 
lines are minima, the dashed line maxima. Up to \QH\/A 2 = 
0.084 the smaller Md yn /A correponds to the deeper minimum 
of the potential, from this value on the larger solution becomes 
the stable configuration. 

fig. where one sees that the r.h.s. of the gap equation 
with if 7^ approaches the if = curve from above. 

To discuss the physical content of the phenomenon just 
described one should fix the characteristic scale A. We 
assume that this value is determined by the problem un- 
der consideration. Its choice can also be motivated by 
the number of Landau levels to be considered. 

In the region m 2 /A 2 <C 1 the four-fermion interaction 
dominates the behaviour of the system. Since their cou- 
pling strength is small, G < G cr , the massless fermions 
behave like almost free particles moving in a weak ex- 
ternal magnetic field, with access to a large number of 
Landau levels, £ 3> 1. This field catalyzes the process of 
fermion-antifcrmion pairing on the energy surface Eq = 
of the LLL. The first minimum localized at to 2 /A 2 <C 1 
is exactly formed by such a (1 + l)-dimensional conden- 
sate. If six- and eight-fermion forces would not act on 
the system, this ground state would be stable: our for- 
mulae as well as the result of paper [j| show clearly that 
a slow increase of the strength if does not wash out 
the condensate from the energy surface Eq = of the 
LLL. However, when the six- and eight-fermion inter- 
actions are present, a slowly increasing magnetic field 
destroys finally this ground state. The new condensate 
has a (3 + l)-dimensional structure similar in every re- 
spect to the standard NJL case with broken chiral sym- 
metry at if = 0, i.e., when the condensate spreads over 
many single fermion states. This is because the increas- 
ing magnetic field enlarges the dynamical fermion mass, 
and scales of order m/A ~ 1 become relevant. At these 
scales the 't Hooft and eight-quark interactions push the 
system to a new regime, where the fermions are not any- 
more free-like particles: they interact strongly with each 
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other and this interaction changes the fermionic spectrum 
and the structure of the ground state in an essential way, 
with all the above mentioned consequences. 

Thus we have obtained not only a correct descrip- 
tion of the well-known physics related with the LLL, but 
have found also a clear signature for the possibly impor- 
tant role played by 't Hooft and eight-quark interactions. 
Namely, in the presence of these interactions the mag- 
netic field can change the condensation zone from the 
zero-energy surface of the LLL to a wide region, spread 
over many Landau levels and vice versa. One can expect 
that hard gamma emissions accompany this process. 

What is the characteristic scale of the magnetic fields 
which can induce such a transition? We obtained the 
value H = 7.3 • 10 13 A 2 Gauss/MeV 2 which actually de- 
pends on the cutoff involved in the problem. For in- 
stance, in hadronic matter it is probably safe to assume 
that A ~ 800 MeV, leading to H = 4.7 • 10 19 G. 



One can indicate several potentially interesting areas 
where this effect may find applications. One of them is 
connected with the recent studies of compact stellar ob- 
jects in presence of strong magnetic fields, in particular 
the young neutron stars, magnetars [25|. The surface 
magnetic fields are observed to be > 10 15 G, but actually 
they can be even much higher in the core region. The 
other area is connected with the electroweak phase tran- 
sition in the early Universe [26], where the strength of 
magnetic fields can reach H ~ 10 24 G. 
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